Probing magnetic order in LLMPO4, M = Ni, Co, Fe and lithium diffusion in Li T FeP04 
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Muon spin relaxation measurements are reported on three members of the Li A MP04 series. The magnetic 
properties of stoichiometric samples with M = Ni, Co, Fe, were investigated at low-temperature. In LiNiPC>4 
we observe different forms of the muon decay asymmetry in the commensurate and incommensurate antiferro- 
magnetic phases, accompanied by a change in the temperature dependence of the muon oscillation frequency. 
In LiCoP04 the form of the muon decay asymmetry indicates that the correlation between layers decreases as 
the Neel temperature is approached from below. LiFePC>4 shows more conventional behaviour, typical for an 
antiferromagnet. Measurements on Li A FeP04 with x = 0.8, 0.9 & 1 show evidence for lithium diffusion below 
~ 250 K and muon diffusion dominating the form of the relaxation at higher temperature. The thermally acti- 
vated form of the observed hopping rate suggests an activation barrier for lithium diffusion of ~ 100 meV and 
a diffusion constant of £)tj ~ 10 -10 — 10~ 9 cm 2 s~' at room temperature. 

PACS numbers: 76.75.+i, 75.50.Ee, 82.47.Aa, 82.56.Lz 



I. INTRODUCTION 

The series of phosphates L1MPO4 crystallize in the or- 
thorhombic olivine structure, with layers of magnetic tran- 
sition metal M ions that are relatively well coupled, mean- 
ing they are intermediate between two- and three-dimensional 
magnetism. 12 The choice of M ion allows the single ion in- 
teractions to be tuned discretely and a range of magnetic be- 
haviour results. In LiNiPC>4 the low-temperature commen- 
surate antiferromagnetic state becomes incommensurate just 
below the bulk ordering temperature. 2 Both LiCoPCU and 
LiNiP04 exhibit magnetoelectric behaviour 3 and the result- 
ing toroidal domains in LiCoP04 have been observed using 
optical measurements^ 

A separate interest in this series of phosphates comes from 
their application as battery cathode materials. This is partic- 
ularly relevant for LiFePC>4 which has a slightly lower cell 
voltage and energy density than the widely used LiCo0 2 , but 
a significantly better lifetime, resistance to thermal runaway, 
and a smaller environmental impact. 5-9 The use of Li A FeP04 
as a battery cathode material leads to questions concerning 
its electrochemical properties and the kinetics of lithium dif- 
fusion, both of which have received considerable study. 9 Both 
calculations and experiment have addressed the activation bar- 
riers for lithium ion and electron conduction, as well as the 
lithium ion diffusion rate4£~— While the calculated values 
have converged there are considerable variations in the results 
of experiments carried out using different techniques. 

In this paper we present a muon-spin relaxation (/jSR) 
investigation of the low-temperature magnetic properties of 
LiMP04 (M = Ni, Co, Fe) and the high-temperature diffusive 
properties of Li A FeP04 (x — 0.8, 0.9, 1.0). As well as being 
a sensitive probe of magnetic ordering, 24,25 /jSR provides a 
means of investigating diffusion processes of both the muon 26 
and other species that perturb its environment. 27 Lithium dif- 
fusion is a process that provides such a perturbation and /jSR 
has now been applied to studying it in a wide range of bat- 
tery cathode materials: Li A [Mni 96Lio,o4]C>4 j 28 ' 29 Lio eTi0 2 , 30 
Li 3 _ x _ v NLN,2J., Li^Co02^^ and LiNiQ 2 ^ Of these, the 



studies on Li A Co0 2 have been the most extensive and found 
that the lithium diffusion rate in this compound is well suited 
to the timescale probed by /jSR. Similar studies on a different 
timescale can be carried out using NMR . 23 - 34 ' 35 

We describe the preparation of the samples and the gen- 
eral /jSR technique in section [Tj] including the details of how 
the low-temperature data were analysed. In sections IIII Al 
IIII Bl and IIII CI we present our /jSR data and analysis. In 
section [IV] we discuss the existing literature on lithium diffu- 
sion in Li A FeP04, describe our higher-temperature /jSR exper- 
iments, data analysis, and results. Our conclusions are sum- 
marized in sectionlVl 



II. EXPERIMENTAL METHOD 

Powders of stoichiometric and Li-deficient Li x MP04 
(M =Ni, Co and Fe; x = 0.8 and 0.9) were synthesized by the 
solid state reaction technique. Starting materials of high purity 
(> 99.99 %) Li 3 P0 4 , NiO, C03O4, Fe 2 3 , and NH 4 H 2 P0 4 
were mixed and sintered in three stages; 175° C for 10 h, 
225° C for 5 h, and 725° C for 24 h. After grinding, they were 
sintered again at 750° C for 24 h. Finally the powders were 
made into rods and sintered at 775° C for 12 h. Single crystals 
of stoichiometric L1MPO4 were grown in a four-mirror optical 
floating-zone furnace (Crystal System Inc.). The growth was 
carried out at a speed of 2-3 mm/h with the feed and seed rods 
counter rotating at 25 rpm. Crystals were grown in an argon 
pressure of \-A atmosphere. 

Spin-polarized positive muons were implanted into the 
samples where they stop at interstitial sites with large elec- 
tronegativity, and decay with a mean lifetime of 2.2 /js. While 
the muons remain within the sample their spin direction is af- 
fected by the local magnetic field at their stopping site, with 
the muon's gyromagnetic ratio = 2k x 135.5 MHz T 1 
being intermediate between those of the electron and pro- 
ton. The muon spin polarization is followed as a function 
of time by measuring the asymmetry in the count rate of de- 
cay positrons, A(t), in two detectors on opposite sides of the 
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FIG. 1: (Color online) Muon decay asymmetry data for L1MPO4: (a) 
M =Ni, (b) M =Co, and (c) M =Fe. Fits to the data are described in 
the text. The initial asymmetries are reduced from the values due to 
energy selection of decay positrons by including data from spatially 
segmented detectors giving better counting statistics. 

sample. 24 Our /jSR experiments were carried out at the Paul 
Scherrer Institute using the General Purpose Surface-muon in- 
strument (GPS) for the low-temperature measurements of the 
stoichiometric crystalline samples and at the ISIS Facility us- 
ing the MuSR spectrometer for the high-temperature measure- 
ments of the Li v FeP04 samples. 

The data shown in Fig.[T]were analysed using the WiMDA 
program^ It was found that similar fitting functions were 
suitable for describing the data on each of the three stoichio- 
metric samples, based on the general form: 

A{t) = A ± e- l 'cos{2% ft +ty)+A ll e- At +A bg . (1) 

The A± term describes an exponentially damped oscillation 
due to a quasistatic magnetic field perpendicular to the muon 
spin direction. The phase could be fixed to zero in both the 
low-temperature commensurate phase of LiNiP04, and below 
Tn in LiFeP04. In LiCoPCU, was found to depend on tem- 
perature. The A || term describes the exponential relaxation for 
muon spins with their direction along that of the local field 
at their stopping site, which are depolarized by spin fluctua- 
tions. The final term describes the temperature-independent 
contribution to the asymmetry from muons stopping outside 



the sample. Just above the magnetic ordering transition there 
is no oscillatory signal and, as is generally the case in para- 
magnets, the data are well described by an exponential relax- 
ation, with rate A. This form of the data is not compatible with 
short-ranged static magnetic order persisting above the long- 
range magnetic ordering transition on the timescales probed 
by muons. We discuss the form of the muon depolarization at 
higher temperatures in section HVl 

In the incommensurate antiferromagnetic phase of 
LiNiPCU we found that Eq. Q] did not provide a satisfactory 
description of the data, even allowing to change from the 
value of zero found to describe the commensurate phase. For 
incommensurate magnetic phases, where the muons sample 
a significant range of magnetic fields, the oscillatory part 
of the muon relaxation function takes the form of a Bessel 
function, 25 so that A(t) can be written as: 

A(f) = A ± e- X '%(2nft) +A\\e- Al +A bg . (2) 

For well-defined oscillation frequencies that varied contin- 
uously below Tn we fitted the temperature dependence to the 
phenomenological function: 

/(r)=/(0)[i-(r/r N )«f , (3) 

where a describes the T —> trend and p describes the trend 
approaching 7n. 

III. LOW-TEMPERATURE RESULTS 

A. LiNiP0 4 

On cooling LiNiP04 orders incommensurately at = 
21.8 K, then orders commensurately at 7ic-c = 20.7 Ki 2 ^ 
Neutron diffraction has characterized both the long-range or- 
dered phases and found diffuse scattering well above T^. 
The trend in the magnetic order parameter approaching the 
breakup of collinear order was found to follow the depen- 
dence expected for a 2D Ising model; 38 an anomalous correla- 
tion between the spin wave spectrum and the incommensurate 
magnetic order 39 and coexisting short- and long-range incom- 
mensurate magnetic order were reported for the intermediate 
phase i 2 - 

In LiNiP04 we can divide the data into three distinct tem- 
perature regions. At the lowest temperatures, below 20.7 K, 
we observe underdamped oscillations at a single frequency, 
with the data shown in Fig.Q](a) well described by EquationQ] 
with = 0. This is consistent with the commensurate mag- 
netic structure determined by neutron diffraction. 2 The tem- 
perature dependence of the oscillation frequency, /, is shown 
in Fig. |2] (a), the linewidth X in (b), and the relaxation rate A 
in (c). The linewidth is relatively small compared to the oscil- 
lation frequency, as is evident from the persistence of the os- 
cillations in the low-temperature data, and only grows slightly 
approaching the commensurate-incommensurate phase transi- 
tion. 

Between 20.7 and 21.8 K the oscillations persist, but their 
form changes from the cosinusoidal form described by equa- 
tion [T] with = to the Bessel function form described by 
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FIG. 2: (Color online) Parameters derived from fitting equations [T] 
and[2]to the juSR data for LiNiPC>4 shown in Fig.[7](a): (a) Oscillation 
frequencies / and /'. (Inset) Reduced oscillation frequency vs re- 
duced temperature (see text) showing the kink at the commensurate- 
incommensurate phase transition, with guides to the eye plotted for 
each of the two phases, (b) Linewidths X and X' . (c) Relaxation rate 
A. 

equation [2] as expected from the incommensurate behaviour 
determined by neutron diffraction. 2 While it is possible to ob- 
tain convergent fits to Eq. Q] with as a free parameter the 
quality of fit is markedly poorer than for Eq.[2]and the differ- 
ence is obvious even to the eye. The oscillation frequency, /', 
linewidth, A/, and relaxation rate A are shown in Fig. [2] (a), 
(b), and (c) respectively. The linewidth grows more rapidly in 
this temperature region which suggests that it is dominated by 
critical fluctuations approaching 7n. 

Fitting the oscillation frequencies shown in Fig.[2](a) using 
Eq. [3j extended to two phases with different p values in each 
phase but a continuous order parameter, leads to the parame- 
ters: T N = 21.76(1) K, /(0) = 26.57(1) MHz, a = 4.22(5), 
Pc = 0.220(3), and Pic = 0.40(5), the last two parameters 
being in the commensurate (C) and incommensurate (IC) or- 
dered phases respectively. The relaxation rate A does not dis- 
play any clear trend in the lowest temperature phase but there 
is a distinct rise in the incommensurate phase. Above the 
data take the exponential form typical of paramagnets where 
the electronic spin fluctuations are fast compared to the distri- 
bution of local fields. The relaxation rate also appears to vary 



FIG. 3: (Color online) Parameters derived from fitting Eq. [TJto the 
/jSR data for LiCoP04 shown in Fig.Q](c). ( a ) Oscillation frequency 
/ with the solid line showing a fit to Eq.[3]and the dashed line show- 
ing the model derived in Ref. 40. (Inset) Phase of oscillating com- 
ponent plotted against the oscillation frequency /. This shows that 
they do not have a monotonic relation, excluding the possibility that 
a time offset is causing the phase offset, (b) Linewidth X. (c) Relax- 
ation rate A. 



critically above T$. 

In the inset to Fig. [2] (a) we show the reduced oscilla- 
tion frequencies plotted against reduced temperature, taking 
f(T = 0) and T N as the fixed points for F = f{T)/f(0) and 
t = (7n — T)/Tn respectively. Trends following the expected 
P values for the two phases are sketched as guides to the eye, 
without the a parameter used in fitting the trend shown in the 
main panel. A kink in the plot of F against t is evident at the 
commensurate-incommensurate phase transition. 

B. LiCoP0 4 

LiCoPCU has been found to be a model example of a mag- 
net intermediate between 2D and 3D exchange coupling. 40 
This gives rise to some unusual critical behaviour, with a neu- 
tron diffraction study reporting that the temperature depen- 
dence of the staggered magnetization follows the form ex- 
pected for the 2D Ising model and that the critical scattering 
above = 22 K follows a 3D Ising form4^ The spin waves 
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have also been studied in detail 41 , showing that as well as the 
dispersions predicted by linear spin wave theory, there is an 
anomalous dispersionless excitation at ~ 1.2 meV that was 
suggested to be related to the magnetoelectric effect in this 
material. 

The form of the raw data for LiCoP04, shown in Fig.[T](b), 
is similar to that in the commensurate phase of LiNiPO/t, and 
can be fitted successfully using Eq. Q] However, it was not 
possible to constrain = all the way up to T^. In such a 
situation, caution is required in distinguishing between a sys- 
tematic shift in the time offset of the raw data (arising from 
an error in determining when muons enter the sample) and a 
phase offset due to the magnetic field distribution in the sam- 
ple. Both possibilities were considered in the data analysis. 
Because the muons were implanted into the sample with their 
initial spin direction rotated relative to the symmetry axes of 
the detector system it was possible to use the geometric phase 
offsets to show that only the phase of the signal was vary- 
ing with temperature. The data were subsequently analyzed 
with (j) as a free parameter, leading to the parameters shown in 
Fig.0 

The oscillation frequency shows a smooth, monotonic 
decrease with increasing temperature and can be fitted to 
the empirical form of Eq. [3] with the parameters: /(0) = 
41.20(1) MHz, T N = 21.72(1) K, a = 4.91(4), and (3 = 
0.299(3). This value of p is smaller than that expected for 
the 3D Ising model (p = 0.326), but considerably larger than 
that for the 2D Ising model (P = 0. 125). Vaknin et al~- com- 
pared the results of their neutron diffraction measurements 
to the analytic Onsager/Yang form 42 M f (r) = M f (0)[l - 
sinh _4 (272D/7 , )] 1/ ^ 8 for the temperature dependence of the 
sublattice magnetization in the 2D Ising model, including an 
additional multiplicative term exp[AE-G/ (T — Tjsj)] to describe 
the interlayer fluctuations relevant near the crossover to 3D 
behaviour. Neither the purely 2D model nor the extended ver- 
sion were able to describe f(T) above 15 K, although the ad- 
ditional multiplicative term does bring the predicted order pa- 
rameter closer to the trend we observe. We plot the coupled 
layer model of Ref.l40las a dashed line in Fig.[3](a) to illustrate 
this difference. Given the effectiveness of the analytic model 
below 15 K and the unusually large value of a (dominated 
by these low-temperature data points), it seems that there is 
a gradual crossover in the effective dimensionality of the sys- 
tem around 15 K that is reflected in the form of the muon 
data. This is consistent with the 3D fluctuations occurring on 
the longer timescale probed by muons at a lower temperature 
than they affect neutron diffraction measurements. 

The linewidth X shown in Fig.[3](b) is slightly larger than 
in LiNiP04 at low temperature and appears to grow in two 
stages as is approached: around 15 K there is a small step 
and above 20 K there is a sharper rise associated with critical 
fluctuations. The relaxation rate A shows almost the opposite 
temperature dependence to A, and is considerably larger than 
in either of the other two samples. There is no sign of a criti- 
cal divergence approaching Tn from below and these observa- 
tions suggest that A is dominated by a quasistatic distribution 
of magnetic fields. Above 7n, A shows a more conventional 
critical divergence. 
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FIG. 4: (Color online) Parameters derived from fitting Eq. [TJto the 
/jSR data for LiFePC>4 shown in Fig. [2(b). (a) Oscillation frequency 
/. (b) Linewidth X. (c) Relaxation rate A. 



The temperature dependence of the phase is plotted in the 
inset of Fig. [3] (a). Below 15 K there is only a slow change 
in but it rapidly changes from — 10° to ~ —80° above 15 K. 
This suggests that, approaching 7n, a magnetic inequivalency 
develops between muon stopping sites in a manner akin to the 
intermediate phase in LiNiPC»4, albeit less pronounced. While 
a weakly incommensurate structure could generate such an ef- 
fect, the sharp increase in occurring at the same temperature 
as both the departure from the quasi-2D trend in f(T) and the 
growth in the linewidth suggests that in LiCoPCU the phase 
shift comes from increasing disorder in the stacking of the 
magnetic layers while approaching from below. 



C. LiFePC-4 

LiFePC>4 shows a commensurate collinear antiferromag- 
netic structure below = 50 K, with neutron diffraction mea- 
surements finding p = 0.27(3). 43 ' 44 The Fe moments are ori- 
entated slightly away from the fo-axis, expected on the basis of 
the crystal symmetry. No evidence for short-range order has 
been observed above 7n, in contrast to the other members of 
this series. 

The parameters derived from fits of equation [T] to the 
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FIG. 5: (Color online) Raw ,uSR data for Li 8 FePO 4 at 140 K with 
fits to Eq.[4]as described in the text. 



LiFePC>4 data are shown in Fig. [4] The temperature depen- 
dence of the frequency and linewidth behave conventionally, 
with f(T) being well described by equation [3] with parame- 
ters: /(0) =45.31(2) MHz, T N = 50.87(7) K, a = 3.66(3), 
and p = 0.381(5). This value of p is close to that expected 
for 3D Heisenberg critical behaviour, as opposed to the value 
of p = 0.27(3) previously estimated from neutron scatter- 
ing measurements . 43 ' 44 The linewidth for T <C 7n is consid- 
erably smaller than the oscillation frequency and it diverges 
smoothly approaching Ts. 

The relaxation rate A behaves differently in LiFePC>4, com- 
pared with the other two samples. There is only a small 
increase around but the primary feature occurs at 30 K, 
well below Associated with this feature is an increase 
in the relaxing amplitude of the signal as the temperature is 
increased. Examining the low-temperature data more care- 
fully allows us to identify a further oscillating component 
with an amplitude around 10 % of the primary oscillating 
component, with an oscillation frequency of ~ 120 MHz. 
This component disappears above 30 K and appears to be re- 
placed by the strongly relaxing term that causes an increase 
in A. We attribute this behavior to a ~ 10 % impurity phase 
that is most likely to be FeP0 4 , since Li v FePC>4 is known to 
form an xLiFeP0 4 :(l — x)FePC>4 binary phase mixture^ and 
r N (FeP0 4 )~ 25 K4£ 
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FIG. 6: (Color online) Parameters derived from fitting Eq. [4] to the 
/jSR data for Li A FeP0 4 . (a) Field distribution width A. (b) Fluctua- 
tion rate v. (Inset) Relaxation rate A for Lig 9FeP0 4 . (c) Fluctuation 
rate v plotted against the inverse temperature to illustrate the ther- 
mally activated behavior below about 250 K. The lines plotted are 
fits to the data in the activated region with parameters described in 
the text. 



IV. HIGH-TEMPERATURE RESULTS FOR Li A FeP0 4 

To investigate the lithium diffusion behavior in Li x FeP04 
we measured three samples with x = 0.8, 0.9, and 1 .0 at tem- 
peratures between 75 and 400 K and at fields of and 0.5 mT. 
By measuring at more than one magnetic field at each tem- 
perature it is possible to get a more reliable determination of 
the fitted parameters, since we have more information on how 
the field distribution experienced by the muon is decoupled by 
the field applied parallel to the initial muon spin polarization. 
Examples of the raw data at the two magnetic fields used are 
shown in Fig. [5] with the fits described below. 

For our high-temperature measurements on Li T FeP04 we 
assume a Gaussian distribution of random local fields due to 
the various magnetic moments present in the sample. For a 



static magnetic system this would lead to a muon depolariza- 
tion described by the Gaussian Kubo-Toyabe function. 4 — Any 
fluctuations present within the muon time window, which may 
be caused by either lithium or muon diffusion, can be treated 
using the strong collision approximation, leading to a dynamic 
Kubo-Toyabe function. 46 Analysis of the data measured at a 
series of fields and temperatures using such a dynamic Kubo- 
Toyabe function proved to be unsuccessful. In studies of 
lithium-containing battery materials it has been usual to mul- 
tiply the dynamic Gaussian Kubo-Toyabe function by an ex- 
ponential relaxation to eliminate any magnetic contribution 
to the relaxationi 27 ' 28 i 31 - 33 This does not lead to reliable fits 
to our raw data either. A consistently better quality of fit 
was obtained by applying Keren's analytic generalization of 
the Abragam function appropriate for /jSR, 47 multiplied by a 
temperature-independent relaxation rate fixed for each sam- 
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pie: 

P z (t) =exp[-r(A,v,C0i,?>]exp(-A/). (4) 

where r(A,v,COi,f) is defined in Ref. 47. The parame- 
ter A describes the quasistatic distribution of field at the 
muon stopping site, v is the temperature-dependent fluctua- 
tion rate, col = y^B is the muon's Larmor precession frequency 
in the applied magnetic field, and A, is due to temperature- 
independent fluctuations. (After initial unconstrained fits had 
been made, A. values were fixed at 0.05, 0.02, and 0.1 MHz 
for the x = 0.8, 0.9, and 1.0 samples respectively.) In the 
x = 0.9 sample we found a strong temperature independent 
relaxation coming from a minority phase which could be sub- 
tracted from the data analysis using A m exp(— At), with the 
values of A shown in the inset to Fig. |6](b). The values of A 
and v obtained from these fits are shown in Fig.|6](a) and (b) 
respectively. 

The A values in Fig. [6] (a) show the trend observed in the 
vast majority of lithium-containing battery materials investi- 
gated to date for x = 0.9, where a low-temperature plateau 
is followed by a smooth decrease to a high-temperature 
plateau i 27 i 28 - 31 i 33 In the x = 0.8 and 1 samples there is a peak 
at around the temperature where the low-temperature plateau 
ends in the x = 0.9 sample. The values of A ~ 0.2 MHz 
are broadly similar to those in Li x CoC>2 27 , LiM^CU^ and 
Lii_ x Nii +x 02*^but smaller than in Li3_ A _ v Nij;Ni^ 

The temperature dependence of v follows a similar trend 
in each sample, with a slight fall from the lowest measured 
temperature to around 100 K, followed by a smooth rise to- 
wards ~ 250 K, and then a sharp drop-off to either the low- 
temperature value, or in x = 0.8, to the value at the peak. 
It seems likely that the change observed below 100 K is 
due to the buildup of magnetic correlations that are not well 
described by our temperature-independent X value. Above 
100 K the thermally activated growth in the hopping rate mir- 
rors that observed in Li x CoC>2 27 and Lii_ T Nii +j: 02,— albeit 
with a different temperature scale. The behavior above the 
peak at ~ 250 K may be related to the onset of muon hopping, 
but this may not be a unique explanation. 

Comparing our results to those obtained previously nat- 
urally leads to the question of whether the phenomena we 
observe are associated with the diffusion of lithium and/or 
muons. The similarity of the temperature dependences of 
both A and v to previous results on other materials indeed 
suggest that they are caused by the same phenomenon. This 
leaves the further question of whether we can obtain quantita- 
tive information about the lithium diffusion from our results. 
While we could not use the dynamic Gaussian Kubo-Toyabe 
function multiplied by an exponential used in Refs. [2ft 
and |3J, the modified Keren function we have employed pro- 
vides the same information and a more robust fit of our data 
over the whole measured temperature range. That v follows 
an activated temperature dependence over a similar tempera- 
ture range to that observed in other materials, as is illustrated 
in Fig. 0(c), strongly suggests that up to around 250 K we can 
assign the change in v to lithium diffusion. Above this temper- 
ature it is likely that the muons become mobile and this results 
in either a drop in A or v as the form of the data changes. 



TABLE I: Comparison of reported estimates for Du and E a obtained 
using different techniques (at room temperature unless noted). A 
more detailed list of E a values is given in Ref.lla. 
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Arrhenius fits to v over the thermally activated region al- 
low us to estimate the energy barriers E a for lithium diffu- 
sion, which forx = 0.8, 0.9, and 1.0 are 130(10), 80(10), and 
60(10) meV respectively. Extrapolating the fits to 300 K for 
comparison with other measured values gives us estimates of 
the lithium hopping rate at room temperature of 2 x 10 6 s _1 
(x = 0.8), 0.8 x 10 6 s" 1 (x = 0.9), and 0.5 x 10 6 s" 1 (x= 1.0). 
(The extrapolation to room temperature introduces an error of 
~ 50 % in these values whereas the individual points within 
the measured range have errors around 10 %.) 

Taking the primary hopping pathway to be along the b- 
axis^ we can further estimate the diffusion constant for 
LiFePCU. The distance travelled in each hop will be b/2 and 
this leads to a diffusion constant Du = b 2 v /A. Given these as- 
sumptions we estimate Du = 1.9 x 10~ 9 cm 2 s _1 for the x = 
0.8 sample, and 7.6 x 10~ 10 cm 2 s _1 and 4.6 x 10~ 10 cm 2 s _1 
for the x — 0.9 and x = 1.0 samples respectively. 

We can compare our estimates for the activation barrier 
and diffusion constant to those derived from other techniques, 
which are summarized in Table |U Most theoretical work and 
experiments find E a for lithium diffusion within the range 
600 - 750 meVj 10 ! 18 ! 22 ! 49 ! 50 although there is both theoretical 
and experimental evidence for E a ~ 100 - 300 meV . 16 ' 19 ' 20 ' 51 
Smaller energy barriers have been suggested for the elec- 
tronic conduction via the hopping of small polarons and it has 
been argued that the polaron hopping may be correlated with 
the hopping of lithium ions , i ' The activation energy of 
^100 meV that we observe suggests that the hopping pro- 
cess observed by the muons is unlikely to be associated with a 
barrier as large as 600 meV and this suggests that another, per- 
haps correlated, process facilitates lithium diffusion at lower 
temperatures. 

The disparity between measurements of Du from differ- 
ent techniques is far greater than that seen for E a , with val- 
ues ranging from 10~ 16 to 10~ 7 cm^- 1 ! 10 ' 12 — Theoretical 
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work and local measurements, such as Mossbauer spec- 
troscopy j 10 ! 11 seem to give larger values of Du ~ 10~ 13 - 
10~ 7 cm 2 r' than bulk measurements, which give Du ~ 
10- 16 - 1(T 10 cm 2 s-' (Ref.QJ). Our estimate lies within 
the overlap of these groups. This suggests that there is a 
difference between microscopic and bulk determinations of 
Du which could result from the effect of the LiFePOzt/FePCU 
phase boundary motion or mesoscopic barriers to lithium dif- 
fusion such as the blocking of diffusion channels by Fe^ de- 
fects and grain boundaries, the latter accentuated by the habit 
of crystallites to be platelets with the b-axis as their shortest 
dimension. 



V. CONCLUSION 

We have used /jSR to provide a new window on both 
the magnetic and diffusive properties of this series of 
olivine phosphates. This has shown how the commensurate- 
incommensurate phase transition in LiNiPCU occurs without 
a discontinuity in the internal field at the muon stopping site; 
how the nature of the fluctuations approaching Tn in LiCoPCU 
are more three-dimensional in the muon time window than 
those found to be quasi-two dimensional in neutron scattering 



measurements; and that the ordering of LiFePCU is more con- 
ventional than the other two materials studied, though again 
the three-dimensional fluctuations are more evident in deter- 
mining the behavior of the order parameter approaching 7n. 
Our measurements also provide a new means of investigating 
the process of lithium diffusion in Li x FeP04, finding a dif- 
fusion constant Du ~ 10 -10 — 10~ 9 cm 2 s~' and an energy 
barrier of E a ~ 100 meV. 

Shortly before we submitted this paper Ref. HH was pub- 
lished reporting analogous measurements of LiFeP04. Two 
oscillating components and a fast relaxing component were 
observed to extend up to Tn suggesting that the higher fre- 
quency component may originate in a metastable muon site or 
nearly degenerate muon sites between which the muon hops. 
The high-temperature data were parameterized slightly differ- 
ently but led to nearly identical values of both Du and E a . 

Acknowledgments 

Parts of this work were performed at the Swiss Muon 
Source, Paul Scherrer Institute, Villigen, CH and the ISIS Fa- 
cility, UK. We thank A. Amato for experimental assistance 
and the EPSRC and STFC (UK) for financial support. 



* Electronic address: peter.baker@stfc.ac.uk 

1 S. Geller and J. L. Durand, Acta Cryst. 13, 325 (1960). 

2 D. Vaknin, J. L. Zarestky, J.-P. Rivera, and H. Schmid, Phys. Rev. 
Lett. 92, 207201 (2004). 

3 I. Kornev, M. Bichurin, J.-P. Rivera, S. Gentil, H. Schmid, A. G. 
M. Jansen, and P. Wyder, Phys. Rev. B 62, 12247 (2000). 

4 B. B. Van Aken, J.-P. Rivera, H. Schmid, and M. Fiebig, Nature 
449, 702 (2007). 

5 A. K. Padhi, K. S. Nanjundaswamy, and J. B. Goodenough, J. 
Electrochem. Soc. 144, 1 188 (1997). 

6 A. Yamada, S. C. Chung, and K. Hinokuma, J. Electrochem. Soc. 
148, A224 (2001). 

7 M. S. Whittingham, Y. Song, S. Lutta, P. Y. Zavalij, and 
N. A. Chernova, J. Mater. Chem. 15, 3362 (2005). 

8 T. Ohzuku and R. Brodd, J. Power Sources 174, 449 (2007). 

9 L.-X. Yuan, Z.-H. Wang, W.X. Zhang, X.-L. Hu, J.-T. Chen, Y- 
H. Huang, and J. B. Goodenough, Energy Environ. Sci. 4, 269 
(2011). 

10 B. Ellis, L. K. Perry, D. H. Ryan, and L. F. Nazar, J. Am. Chem. 
Soc. 128, 11416 (2006). 

11 J. L. Dodd, I. Halevy, and B. Fultz, J. Phys. Chem. C 111, 1563 
(2007). 

12 P. P. Prosini, M. Lisi, D. Zane, and M. Pasquali, Solid State Ionics 
148, 45 (2002). 

13 A. V. Churikov, A. V. Ivanishchev, I. A. Ivanishcheva, V. O. Sy- 
cheva, N. R. Khasanova, and E. V. Antipov, Electrochem. Acta 
55, 2939 (2010). 

14 S. Franger, F. Le Cras, C. Bourbon, and H. Rouault, Electrochem. 
Solid-State Lett. 5, A231 (2002). 

15 D. Y. W. Yu, C. Fietzek, W. Weydanz, K. Donoue, T. Inoue, 
H. Kurokawa, and S. Fujitani, J. Electrochem. Soc. 154, A253 
(2007). 

16 D. Morgan, A. Van der Ven, and G. Ceder, Electrochem. Sol. -Stat. 



Lett. 7, A30 (2004). 

17 M. S. Islam, D. J. Driscoll, C. A. J. Fisher, and P. R. Slater, Chem. 
Mater. 2005, 5085. 

18 R. Amin, J. Maier, P. Balaya, D. P. Chen, and C. T. Lin, Solid 
State Ionics 179, 1683 (2008). 

19 M. Takahashi, S. Tobishima, K. Takei, and Y. Sakurai, Solid State 
Ionics 148, 283 (2002). 

20 C. Wang and J. Hong, Electrochem. Solid-State Lett. 10 A65 
(2007). 

21 T. Maxisch, F. Zhou, and G. Ceder, Phys. Rev. B 73, 104301 
(2006). 

22 K. Hoang and M. Johannes. TarXiv: 1 105.3492 (unpublished). 

23 J. Cabana, J. Shirakawa, G. Chen, T. J. Richardson, and C. P. 
Grey, Chem. Mater. 22, 1249 (2010). 

24 S. J. Blundell, Contemp. Phys. 40, 175 (1999). 

25 A. Yaouanc and P. D. de Reotier, Muon Spin Rotation, Relax- 
ation, and Resonance (Oxford University Press, Oxford, 201 1). 

26 R. Kadono, J. Imazato, T. Matsuzaki, K. Nishiyama, K. 
Nagamine, T. Yamazaki, D. Richter, and J.-M. Welter, Phys. Rev. 
B 39, 23 (1989). 

27 J. Sugiyama, K. Mukai, Y. Ikedo, H. Nozaki, M. Mansson, and 
I. Watanabe, Phys. Rev. Lett. 103, 147601 (2009). 

28 C. T. Kaiser, V. W. J. Verhoeven, P. C. M. Gubbens, F. M. Mulder, 
I. de Schepper, A. Yaouanc, P. Dalmas de Reotier, S. P. Cottrell, 
E. M. Kelder, and J. Schoonman, Phys. Rev. B 62, R9236 (2000). 

29 M. J. Ariza, D. J. Jones, J. Roziere, J. S. Lord, and D. Ravot, J. 
Phys. Chem. B, 2003, 6003. 

30 P. C. M. Gubbens, M. Wagemaker, S. Sakarya, M. Blaauw, 
A. Yaouanc, P. Dalmas de Reotier, and S. P. Cottrell, Solid State 
Ionics 177, 145 (2006). 

31 A. S. Powell, J. S. Lord, D. H. Gregory, and J. J. Titman, J. Phys. 
Chem. C 113, 20758 (2009). 

32 K. Mukai, J. Sugiyama, Y. Ikedo, H. Nozaki, K. Shimomura, 



8 



K. Nishiyama, K. Ariyoshi, and T. Ohzuku, J. Power Sources 174, 
711 (2007). 

33 J. Sugiyama, Y. Ikedo, K. Mukai, H. Nozaki, M. Mansson, O. 
Ofer, M. Harada, K. Kamazawa, Y. Miyake, J. H. Brewer, E. J. 
Ansaldo, K. H. Chow, I. Watanabe, and T. Ohzuku, Phys. Rev. B 
82, 224412 (2010). 

34 V. W. J. Verhoeven, I. M. de Schepper, G. Nachtegaal, A. P. M. 
Kentgens, E. M. Kelder, J. Schoonman, and F. M. Mulder, Phys. 
Rev. Lett. 86, 4314(2001). 

35 K. Nakamura, H. Ohno, K. Okamura, Y. Michihiro, T. Moriga, 
I. Nakabayashi, and T. Kanashiro, Solid State Ionics 177, 821 
(2006). 

36 F. L. Pratt, Physica B 289-290, 710 (2000). 

37 Yu. N. Kharchenko, N. F. Kharcheno, M. Baran, and R. Szym- 
czak, Low Temp. Phys. 29, 579 (2003). 

38 D. Vaknin, J. L. Zarestky, J. E. Ostenson, B. C. Chakoumakos, 
A. Goni, P. J. Pagliuso, T. Rojo, and G. E. Barberis, Phys. Rev. B 
60, 1100(1999). 

39 T. B. S. Jensen, N. B. Christensen, M. Kenzelmann, 
H. M. R0nnow, C. Niedermayer, N. H. Andersen, K. Lefmann, 
M. Jimenez-Ruiz, F. Demmel, J. Li, J. L. Zarestky, and D. Vaknin, 
Phys. Rev. B 79, 092413 (2009). 

40 D. Vaknin, J. L. Zarestky, L. L. Miller, J.-P. Rivera, and 
H. Schmid, Phys. Rev. B 65, 224414 (2002). 

41 W. Tian, J. Li, J. W. Lynn, J. L. Zarestky, and D. Vaknin, Phys. 



Rev. B 78, 184429 (2008). 

42 L. Onsager, Phys. Rev. 65, 117 (1944); C. N. Yang, ibid. 85, 809 
(1952). 

43 G. Rousse, J. Rodriguez-Carvajal, S. Patoux, and C. Masquelier, 
Chem. Mater. 2003, 4082. 

44 J. Li, V. O. Garlea, J. L. Zarestky, and D. Vaknin, Phys. Rev. B 73, 
024410 (2006). 

45 P. D. Battle, A. K. Cheetham, C. Gleitzer, W. T. A. Harrison, G. 
J. Long, and G. Longworth, J. Phys. C, 15, L919 (1982). 

46 R. S. Hayano, Y. J. Uemura, J. Imazato, N. Nishida, T. Yamazaki, 
and R. Kubo, Phys. Rev. B 20, 850 (1979). 

47 A. Keren, Phys. Rev. B 50, 10039 (1994). 

48 S. Nishimura, G. Kobayashi, K. Ohoyama, R. Kanno, M. 
Yashima, and A. Yamada, Nature Materials 7, 707 (2008). 

49 C. Delacourt, P. Poizot, J.-M. Tarascon and C. Masquelier, Nature 
Mater. 4, 254 (2005). 

50 J. Molenda, W. Ojczyk, K. Swierczek, W. Zajac, F. Krok, J. Dy- 
gas, and R.-S. Liu, Solid State Ionics 177, 2617 (2006). 

51 S. Q. Shi, L. J. Liu, C. Y. Ouyang, D. S. Wang, Z. X. Wang, L. Q. 
Chen, and X. J. Huang, Phys. Rev. B 68, 195108 (2003). 

52 J. Sugiyama, H. Nozaki, M. Harada, K. Kamazawa, O. Ofer, M. 
Mansson, J. H. Brewer, E. J. Ansaldo, K. H. Chow, Y. Ikedo, Y. 
Miyake, K. Ohishi, I. Watanabe, G. Kobayashi, and R. Kanno, 
Phys. Rev. B 84, 054430 (201 1). 



